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The features of net baryon productions and collective flow in relativistic heavy-ion collisions at 
energies reached at the CERN Large Hadron Collider (LHC), BNL Relativistic Heavy Ion Collider 
(RHIC) , CERN Super Proton Synchrotron (SPS) and BNL Alternating Gradient Synchrotron 
(AGS) with the model of Non-Uniform Flow Model(NUFM) are systematically studied in this 
paper. Especially we predict the feature of net baryon productions and collective flow at LHC 
^/i^T=5500 GeV basing on the detailed study of that at RHIC v /s^77=62.4 and 200GeV. The 
dependencies of the features of baryon stopping and collective flow on the collision energies and 
centralities are investigated. 
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I. INTRODUCTION 

Over the last two decades, relativistic heavy ion colli- 
sions [H, 0| have been studied experimentally at increas- 
ingly higher center-of-mass energies at the Brookhaven 
Alternating Gradient Synchrotron AGS sjvjv < 
5 GeV), the CERN Super Proton Synchrotron 
SPS (^/sjvat < 20 GeV) and the Brookhaven Relativistic 
Heavy Ion Collider RHIC (y/sW^ < 200 GeV). As dis- 
cussed in this article, the data collected in these experi- 
ments display remarkable generic trends as a function of 
system size and kinematic variables. The Large Hadron 
Collider LHC at CERN will study heavy ion collisions 
at a center-of-mass energy y / saT/v = 5.5 TeV, which is 
a factor 27 higher than the maximal collision energy at 
RHIC. This is an even larger increase in center of mass 
energy than the factor 10 in going from the CERN SPS to 
BNL RHIC. It leads a significant extension of the kine- 
matic range in longitudinal rapidity and transverse mo- 
mentum. The collectivity of high energy density matter 
is one of the important properties to understand high- 
energy heavy- ion collisions 3-5]. It is also challenging 
to understand how collectivity is generated during colli- 
sions. 

There has been a lot of work in recent years on thermal 
and collective flow model calculations d-[I3 of heavy ion 
collisions to RHIC data and extrapolating them to the 
higher LHC energies. Here we should mention some kinds 
of models of thermal and collective flow. The first one is 
the spherically-expanding source model that may be ex- 
pected to approximate the fireball of an isotropic thermal 
distribution created in lower-energy collisions. 

As the collision energy increases, stronger longitudinal 
flow is formed which leads to a cylindrical geometry ac- 
cording to the second kind model [lg, iLJJ . It accounted 
for the anisotropy of longitudinal and transverse direc- 
tion by adding the contribution from a set of fireballs 



with centers located uniformly in the rapidity region in 
the longitudinal direction. It can account for the wider 
rapidity distribution at AGS and SPS when comparing 
to the prediction of pure thermal isotropic model. 

Bjorken [l8[ postulated that the rapidity distribu- 
tion of produced particles establishes a plateau at mid- 
rapidity which has been formulated for asymptotically 
high energies. It is well known that collisions at available 
heavy-ion energy regions of AGS , SPS and RHIC are nei- 
ther fully stopped nor fully transparent (h^-Ht} , although 
a significant degree of transparency is observed. But the 
central plateau structure becomes more and more obvi- 
ous as the collision energy increases to SPS and RHIC. 

As the collision energy increases to LHC, which is a 
factor 27 higher than the maximal collision energy at 
RHIC, the kinematic range in the longitudinal direction 
will increase considerably and the net-baryon density will 
decrease quickly at mid-rapidity. It seems reasonable to 
realize that the plateau proposed by Bjorken [l8| at mid- 
rapidity at LHC energy region has been established. For 
the net baryon distributions, Ref. [28| realized the colli- 
sion of high-energy heavy ions can be divided into two 
different energy regions: the baryon-free quark gluon 
plasma( QGP) region (or the pure QGP region) with 
y/s > lOOGeV per nucleon, and the baryon-rich QGP re- 
gion (or the "stopping" region) with y/s ~ 5 — lOGeV 
per nucleon, which corresponds to about many tens of 
GeV per projectile nucleon in the laboratory system. In 
the baryon-free QGP region, we need to know the nuclear 
stopping power to determine whether the beam baryons 
and the target baryons will recede away from the center 
of mass without being completely stopped, leaving behind 
QGP with very little baryon content. 

The NUFM (Non-Uniform Flow Model) [29l - [33 | realized 
that the fireballs keep some memory on the motion of the 
incident nuclei, and therefore the distribution of fireballs, 
instead of being uniform in the longitudinal direction, is 
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more concentrated in the motion direction of the incident 
nuclei, i.e. more dense at large absolute value of rapid- 
ity. It will not only lead to anisotropy in longitudinal- 
transverse directions, but also render the fireballs (es- 
pecially for those baryons) distributing non- uniformly in 
the longitudinal direction. NUFM [29l - l34l ] may analyze the 
central dip of baryon rapidity distribution by assuming 
that the centers of fireballs are distributed non-uniformly 
in the longitudinal phase space. 

This paper is organized as follows. In Sec. 2 we give a 
brief review the Non- Uniform Flow Model in the longitu- 
dinal direction. The comparison and analysis of baryon 
distribution of AGS , SPS , RHIC and LHC with the results 
of the model given in Sec. 3. Section 4 gives a summary 
and conclusions. 



II. NON-UNIFORM FLOW MODEL NUFM 



The NUFM model we considered 29-34] contains three 
distinct assumptions: 

1. It is argued that the transparency /stopping of rela- 
tivistic heavy-ion collisions should be taken into account 
more carefully. A more reasonable assumption is that the 
fireballs keep some memory on the motion of the incident 
nuclei, and therefore the distribution of fireballs, instead 
of being uniform in the longitudinal direction, is more 
concentrated in the direction of motion of the incident 
nuclei, i.e. more dense at large absolute value of rapid- 
ity. It will not only lead to anisotropy in longitudinal- 
transverse directions, but also render the fireballs (es- 
pecially for those baryons) distributing non-uniformly in 
the longitudinal direction. 

2. The freeze-out temperatures are assumed to be 
about the same around 120 MeV whether it is at higher 
LHC or at lower AGS energy region. Since the temper- 
ature at freeze-out exceeds 100 MeV, the Boltzmann ap- 
proximation seems reasonable to study LHC at freeze- 
out. 

3. In order to express the non- uniformity of flow in 
the longitudinal direction, an ellipticity parameter e is 
introduced through a geometrical parametrization. For 
the central collisions, the nuclear stopping can be studied 
by the range of rapidity of emission source in the center- 
of-mass system. 

We have previously used NUFM to study the net pro- 
ton rapidity among AGS, SPS and RHIC energy re- 
gions [20| . But for the RHIC energy regions, we made 
an earlier error [301 ] to predict the distributions of net 
proton distributions since we neglected the effects of the 
baryon number conservation. Therefore, it is necessary 
to reanalyze the features of net proton rapidity distribu- 
tions among AGS to RHIC by taking into account the 
baryon number conservation. It is found that when we 
consider the baryon number conservation, the features of 
the distributions at RHIC are completely different from 
the results given before [3(|, especially at large absolute 
rapidity region. On the other hands, with the run of 



forthcoming LHC, the predictions of the features of net 
proton rapidity distributions at LHC are also important. 
We will restudy the features of net proton rapidity distri- 
butions among AGS to RHIC by using NUFM, and make 
prediction for the features of forthcoming LHC in this 
paper. In the following, we will firstly make a simple 
introduction to the NUFM. 

A parametrization for such a non-uniform distribution 
can be obtained by using an ellipse-like picture on emis- 
sion angle distribution. In this scenario, the emission 
angle is: 



= tan 1 (etan9). 



(1) 



Here, the induced parameter e(0 < e < 1) repre- 
sents the ellipticity of the introduced ellipse which de- 
scribes the non-uniform of fireball distribution in the lon- 
gitudinal distribution. The detailed discussions of the 
NUFM were given by ref . [2!| . The rapidity distribution 
of NUFM is: 
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p(y e )dy e (l + 2T + 2T 2 )e~ 1 / r 1 

(2) 

y e o and e are the important parameters in this paper, y e Q 
is the rapidity limit which confines the rapidity interval 
of longitudinal flow and e can describe the non-uniform 
in the longitudinal direction of the collective flow. In eq. 
(2), 



T/mcosh(y -y e ), 



(3) 



m is the mass of produced particle, T is the temperature 
parameter, and y e is the rapidity of collective flow, and 



P{Ve) = 



1 + sinh (y e ) 
1 + e 2 smh 2 (y e 



(4) 



is the flow distribution function in the longitudinal direc- 
tion, and e is a parameter which represents the ellipticity 
of the introduced ellipse describing the non-uniform of 
fireball distribution in the longitudinal direction. It may 
be figured out from Eq.2 to Eq.4 that the larger the pa- 
rameter e, the flatter the distribution function p(y e ), the 
more uniform the longitudinal flow distribution. When 
e =>■ 1, the longitudinal flow distribution is completely 
uniform p(y e ) 1 and returns to uniform flow. The 
other important parameter y e o describes the kinematic 
region and can determine the width of the distribution. 

In order to discuss the dependence of of velocity of col- 
lective flow on collision energy in the central mass(CM) 
system, we give a calculation of the average velocity in 
the longitudinal direction as < /?£ >= tanh(y e o/2) and 
< /?7 >l, where 7 = \j \]\— < f3 L > 2 is the lorentz fac- 
tor. Therefore y e Q can also determine the average velocity 
of collective flow in the longitudinal direction. 
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III. THE NET PROTON DISTRIBUTIONS AT 
THE WHOLE AGS TO LHC ENERGY REGIONS 

We use the form of the NUFM model as described in 
Section II and fit the experimental data with the param- 
eters y e o and e that have been assumed to be different 
for different energies, as given in Table I. The systemat- 
ics of these parameters provide useful information on the 
collective flow of baryons in these reactions. 

Comparing with NUFM calculation before [3(| , we con- 
sider the influence of baryon number conservation at this 
time when discussing the distributions. Fig. 1(a) shows 
net-proton rapidity distributions measured at AGS , and 
SPS energies. Fig. 1(b) shows net-proton rapidity dis- 
tributions of the top 5% central collisions measured at 
RHIC y/s^ = 62.4 and ^/sWn = 200 GeV , respec- 
tively. The solid lines are our NUFM calculation results 
from AGS to RHIC and the dotted line is the calcu- 
lation result for that of LHC . It can be seen from the 
Fig.l that NUFM model can fit the experimental results 
from AGS to RHIC , and reproduce central dip of the 
rapidity distribution of the proton at SPS and RHIC in 
agreement with the experimental findings. y e o is approx- 
imately equal to the half width of fit distribution. In the 
sense, the parameter y e o represent the kinetic region of 
collective flow in the longitudinal direction. The param- 
eter of T is chosen to be 0.12 GeV. 

The features of non-uniform flow distributions show 
strong energy dependence from AGS to RHIC . For ex- 
ample: at AGS (e = 0.82, ^i a b =10.8GeV), the net pro- 
ton distribution has a peak at mid-rapidity, and the dis- 
tribution is narrower than that of the other two ener- 
gies. The collective flow is approximately uniform. While 
at SPS (e = 0.61) a dip begins to show in the mid- 
dle of rapidity distribution. While 
62.4GeV and 
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at RHIC snn = 
200GeV the distributions show 



deep dip and the non-uniform parameter e take 0.34 and 
0.31, respectively. According to our calculation as the col- 
lision energy increases, the net baryon distributions be- 
come wide for the whole rapidity distribution and the 
net baryon densities become small at the middle rapidity 
(y ~ 0) region. 

We also speculate the feature of non-uniform flow dis- 
tributions at LHC in fig. 1(b). As mentioned before, 
although collisions at available heavy-ion energy regions 
of AGS, SPS and RHIC are neither fully stopped nor fully 
transparent, but the plateau structure becomes more and 
more obvious as the collision energy increases to SPS and 
RHIC. It leads a significant extension of the kinematic 
range (y e o)in longitudinal rapidity and the net-baryon 
distribution at the central rapidity region decreases at 
LHC. It seems reasonable that we conjecture the kine- 
matic range {y e o) at LHC ^snn = 5500 GeV approaches 
to the incident beam rapidity y p , which is about 1.5 times 
of that at RHIC ^/sWn = 200 GeV (y e0 = 5.32). 

By making a analysis of the dependence of dN/dy \ v —q 
on incident beam rapidity y p from SPS (^/snn = 
17.2GeV) to RHIC {y/sWN=62A and 200 GeV) experi- 
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FIG. 1: The net proton distribution at AGS and 
SPS are shown in Fig. 1(a). the net proton distribution at 
RHIC yfsNN = 62.4GeV and ^s^n = 200GeV are shown in 
Fig. 1(b). The experimental results come from [13^31, the 
dotted line which is predicted by NUFM for LHC is shown 
in Fig. 1(b). 



ments, in which the rapidity distribution obviously show 
central dip feature, we can provide a relationship between 
dN/dy \ y =o and incident beam rapidity y p 



dN/dy | y=0 = 51.0 - 22.0 • \og(y p ) 



(5) 



shown in Fig. 2. According to the speculation, the mag- 
nitude of rapidity density at central rapidity y at 
LHC is about 3.63 which is about 1/4 times of that at 
RHIC 200GeV. We can get e = 0.19 to fit the LHC dis- 
tribution by using the NUFM and the baryon number 
conservation law. 

At LHC , a broad dip in the middle of rapidity region 
has developed spanning several units of rapidity, indicat- 
ing that collisions are quite transparent at LHC energy 
region. According to our study, e = 0.19 at LHC gives 
a more obvious non-uniform feature than that of AGS , 
SPS and RHIC energy region, and the detailed results 
are shown in Table. 1. 

From Fig. 3, we know that as the incident energy in- 
creases, the longitudinal flow distribution becomes more 
non-uniform, e = 0.19 at LHC v /snn = 5500GeV is 
smaller than e = 0.31 at RHIC snn = 200GeV, and 
e = 0.82 at AGS (E lah =10.8GeV). The central rapidity 
density at AGS is the largest in the whole AGS , SPS , 
RHIC and LHC energy regions in Fig. 3. 

Figure 4 shows proton rapidity distribution at differ- 
ent collision energies at AGS and SPS, and the solid lines 
are the calculation results. From Fig. 4 we know that pro- 
ton distribution shows uniform distribution feature in the 
longitudinal direction at the AGS(2Gev-8GeV) according 
to NUFM. But for Pb + Pb interactions (158GeV) at SPS, 




FIG. 2: The dependence of central rapidity density 
dN/dy \ y= o on incident beam rapidity y p from SPS (^/snn = 
17.2GeV) to RHIC (^s^iv=62.4 and 200 GeV) experiments. 
The solid circles are from experimental results, the solid tri- 
angle is the speculating . The real line is the fit curve of 
Eq.5. 



Figure 5 shows proton distribution at different colli- 
sion systems AGS and SPS. It is found that from Fig. 
5 that e = 1 and y e o— 1.411 for heavier collision system 
(Au + Au), but e = 0.72 and y e0 =1.609 for lighter col- 
lision system (Si + Al) at AGS. It is suggested that the 
lighter the collision system, the more non-uniform the 
distribution in the longitudinal direction, the larger the 
kinematical limitation. The same situation is shown at 
SPS comparing with AGS. 

From the calculation, we find that y e o determines the 
width of distribution and confines the flow kinetics re- 
gions. < /?7 > L 

It is found that the depth of the central dip of the net 
baryon distributions depends on the magnitude of the 
parameter e that describes the non-uniformity of longi- 
tudinal flow. 

The stopping may be estimated from the rapidity loss 
experienced by the baryons in the colliding nuclei. If 
incoming beam baryons have rapidity y p relative to the 
CM, the average rapidity loss of net proton is 

< Sy >= y p -<y> (6) 
where < y > is the average rapidity of net proton. 
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FIG. 3: The flow distribution function of net proton in the 
longitudinal direction in the whole AGS , SPS , RHIC and 
LHC energy regions 



e = 0.61 shows non-uniform distribution feature in the 
longitudinal direction. 



< V >= 



part JO 



Up 



ydy 



dN 



B-B 



(y) 



dy 



(7) 



where N valt is participant nucleon number. y p is rapidity 
of incoming beam baryons relative to the CM . The < y > 
is given by 



dn 



< y >= 



Io Vp ydy 
Io Vp dy% 

where dn/dy is given by NUFM . 



(8) 



TABLE I: The different parameters of net proton distribution 
by using NUFM from AGS to LHC 



E ub or (GeV) 


y P 


<S y > 


</?7 >L 


e 


yeo 


Ei ab = 2 (Au+Au AGS) 


0.6951 


0.3519 


0.3255 


1.0 


0.648 


Ei ab = 4 (Au+Au AGS) 


1.0647 


0.5391 


0.4653 


1.0 


0.910 


E lab = 6 (Au+Au AGS) 


1.2714 


0.6332 


0.5897 


1.0 


1.124 


E lab = 8 (Au+Au AGS) 


1.4166 


0.6997 


0.6189 


1.0 


1.168 


E tab = 10.8 (Au+Au AGS) 


1.5674 


0.9499 


0.6967 


0.82 


1.300 


E lah = 14.6 (Si+Al AGS) 


1.7186 


0.7989 


0.7256 


0.72 


1.684 


Ei ab = 158 (Pb+Pb SPS) 


2.9112 


1.6774 


1.4558 


0.61 


2.340 


E lab = 200 (S+S SPS) 


3.0283 


1.1336 


1.6542 


0.554 


2.960 


^/s^ = 62.4 (Au+Au RHIC) 


4.197 


1.9528 


3.2682 


0.34 


4.860 


ysiiT = 200 (Au+Au RHIC) 


5.36 


2.3021 


7.7894 


0.31 


5.320 


yfs^ = 5500 (Pb+Pb LHC) 


8.4669 


3.5724 


32.7912 


0.19 


7.880 



From Fig. 6(a) and Table. 1, we know that from AGS to 
SPS, average rapidity loss < Sy > increases linearly with 
y p . When discussing at RHIC , we study the average ra- 
pidity loss at v / snn = 62.4 and 200 GeV, a new linear 
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FIG. 4: Proton rapidity distributions for Au + Au interac- 
tion at AGS. The experimental data is from Ref. 22-27J an d 
the solid lines are the calculation results. The whole fitted 
parameters e and quanta y e o are given by Table.! 
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FIG. 6: The dependence of average rapidity loss < Sy > (a) 
and < /?7 >l (b) on incident proton rapidity in the whole 
AGS , SPS , RHIC and LHC energy regions. 



power at LHC. The dependence of < (3-f >l (b) on inci- 
dent proton rapidity in the whole AGS , SPS , RHIC and 
LHC energy regions are shown in Fig. 6(b). we can find 
a kind of Log increasing dependence of < /3^y >l (b) on 
incident proton rapidity. 



>.60 

73 



40 



20 



RHIC Au+Au 



200 GeV 



-4 




FIG. 5: Proton rapidity distributions for Au + Au (10.8GeV) 
and Si + Al (14.6GeV) interactions at AGS, and for Pb + 
Pb (I58GeV) and S + S (200GeV) interactions at SPS. The 
experimental data is from Ref. [22H27I ] , and the solid lines 
are the calculation results. The whole fitted parameters e 
and quanta y e o are given by Table. 1 



increasing relationship is established from SPS to RHIC , 
but begins to increase slowly and deviates from that of 
from AGS to SPS. We also predict the nuclear stopping 



FIG. 7: Rapidity distribution of net baryons in Au + Au 
collisions at RHIC energy of ^/sjvjv = 200GeV are compared 
with preliminary BRAHMS net baryon data [3^] for different 
centralities of 0%-10%, 10%-20%, 20%-40%, and 40%-60%. 

As shown in Figure 7, NUFM can fit the net baryon dis- 
tribution at different centralities of 0%-10%, 10%-20%, 
20%-40%, and 40%-60% at RHIC = 200 GeV. 

From Figure 8 and Table. 2, we know that as the cen- 
trality increases, the kinematic region and average veloc- 
ity < /?7l > in the longitudinal direction increases and 
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TABLE II: The fit parameters of net proton distribution for 
different centralities of 0%-10%, 10%-20%, 20%-40%, and 
40%-60% by using NUFM at RHIC y/s^ = 200 GeV. 



centrality at snn = 200 GeV 


2/e0 


< Sy > 


< /?7 >L 


e 


0%-10% (Au+Au RHIC) 


5.320 


2.318 


7.113 


0.31 


10%-20% (Au+Au RHIC) 


4.699 


2.575 


5.195 


0.31 


20%-40% (Au+Au RHIC) 


4.239 


2.822 


4.106 


0.31 


40%-60%(Au+Au RHIC) 


4.199 


2.844 


4.022 


0.31 



the distribution becomes wide. On the other hands, the 
stopping power diminishes as the centrality increases. It 
is surprise to find that the non-uniformity e keeps un- 
changed. 
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FIG. 8: The dependence of average rapidity loss < Sy > 
(Fig. 8a), < /?7 >l (Fig8.b) and kinematic region y e o (Fig8.c) 
on the collision centrality at RHIC ^/snn = 200 GeV. 



IV. SUMMARY AND CONCLUSION 

Net proton rapidity distributions have been measured 
by several experiments at different energies from AGS to 
RHIC . The compiled data are shown in Fig. 1. The net 
proton rapidity distributions are reconstructed among 
the AGS , SPS ,RHIC and LHC energy regions by us- 
ing NUFM in this work. We can predict the distribution 
feature in the fragmentation region of the net proton dis- 
tributions at RHIC although RHIC [lj| only provided the 
multiplicity distribution of net protons at the central ra- 
pidity region. While at RHIC ^/snn — 62.4GeV and 
tJsnn — 200GeV the distributions show deep dip and 
the non-uniform parameter e take 0.34 and 0.31, respec- 



tively. According to our calculation as the collision en- 
ergy increases, the net baryon distributions become wide 
for the whole rapidity distribution and the net baryon 
densities become small at the middle rapidity (y pa 0) 
region. 

The features of non-uniform flow distributions show 
strong energy dependence from AGS to RHIC . For ex- 
ample: the net proton distribution at AGS has a peak at 
mid-rapidity, and the distribution is narrower than that 
of the other two energies. The collective flow is approx- 
imately uniform. While at SPS a dip begins to show in 
the middle of rapidity distribution. The distributions at 
RHIC V^vjv = 62.4GeV and ^/s^T/v = 200GeV show the 
non-uniform feature of deep dip. It is found that the dis- 
tributions become wider and wider for the whole rapidity 
distribution and the densities of net baryon in the middle 
of rapidity region (y ps 0) become smaller and smaller. 

Here,we should mention that quite a few theoretical 
models [3rj| - l40| can give equally good representation of 
the data of particle productions in relativistic heavy ion 
collisions. These models give some different physical pic- 
tures for the research. In Ref. [38] in order to determine 
whether a pure quark-gluon plasma with no net baryon 
density could be formed in the central rapidity region 
in relativistic heavy ion collisions, Wong [38| estimated 
the baryon distribution by using a Glauber-type multiple 
collisions in which the nucleons of one nucleus degrade in 
energy as they make collisions with nucleons in the other 
nucleus. It was found that in the head-on collision of 
two heavy nuclei(^4 > 100), the baryon rapidity rapidity 
distributions have broad peaks and extend well. 

The plateau structure becomes more and more obvi- 
ous as the collision energy increases to RHIC although 
collisions at available heavy-ion energy regions of AGS, 
SPS and RHIC are neither fully stopped nor fully trans- 
parent. It leads a significant extension of the kinematic 
range (?/ e o)in longitudinal rapidity and the net-baryon 
distribution at the central rapidity region decreases at 
LHC. We can study the feature of net baryon distribu- 
tions at LHC by using the NUFM and the baryon number 
conservation law. 

Detailed energy dependence of the net baryon distribu- 
tion among AGS , SPS and RHIC shows a clear transition 
from the baryon stopping region to the baryon transpar- 
ent region. It is found that from AGS to SPS , aver- 
age rapidity loss < Sy > increases linearly with y p , but 
begins to increase slowly and deviates from the linear 
relationship when at RHIC and LHC . It may suggests 
the difference of the interaction mechanism in RHIC and 
LHC from AGS and SPS . The detailed study of net pro- 
ton rapidity distributions from AGS to LHC will deepen 
our study of the relativistic heavy-ion collisions. 

It is found that the transparency of relativistic heavy- 
ion collisions increases as collision energy increases, i.e. 
the higher the collision energy, the more transparent of 
the collision system by analyzing the proton rapidity dis- 
tribution. The phase space of heavy collision system 
is nearly completely uniform in the longitudinal direc- 
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tion at AGS . The phase space of proton distributes non- 
uniformly in the longitudinal direction at SPS and RHIC . 
At LHC , a broad dip in the middle of rapidity region has 
developed spanning several units of rapidity, indicating 
that collisions are quite transparent at LHC energy re- 
gion. According to our study, e = 0.19 at LHC gives 
a more obvious non-uniform feature than that of AGS , 
SPS and RHIC energy region. By reanalyzing RHIC , w e 
obtain a wider rapidity distribution than that of Ref [30j ■ 
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